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Using high resolution x-ray scattering, we have measured the structure 



■»^j- ■ of the Qi CDW in Ta-doped NbSe3. Detailed line shape analysis of the 

data demonstrates that two length scales are required to describe the phase- 
phase correlation function. Phase fluctuations with wavelengths less than a 
new length scale a are suppressed and this a is identified with the amplitude 

o" 

O ■ coherence length. We find that £ a * = 34.4 ± 10. 3A. Implications for the 



physical mechanisms responsible for pinning are discussed. 
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During the past twenty years, the influence of random disorder on phase transitions has 
been studied extensively Disorder can be formally described as randomness in either the 
interaction strength [|TJ or the field conjugate to the order parameter 0. It is widely believed 
that a small amount of randomness in the interactions is irrelevant. Random fields, on the 
other hand, have dramatic consequences. In their seminal paper, Imry and Ma suggested 
that a random field should cause the lower marginal dimensionality dt to rise from 2 to 4 for 
systems with continuous symmetry. Here, dt is the dimensionality below which the system 
cannot sustain long-range order (LRO) at finite temperatures |§. This loss of LRO has been 
observed in a wide variety of systems. In particular, the structures of the charge-density 
waves (CDW's) found in the quasi-one- dimensional metal NbSe 3 do not exhibit LRO ||,[J. 

For mathematical simplicity, the quenched random field is frequently assumed to be 
a time-independent random Gaussian variable satisfying (h(r)) = and (h(r)h(r')) = 
nih,Q S d (r — r') where n« is the number density of impurity atoms and ho gives the strength of 
the pinning interaction. In the general case, scattering from a random structure is charac- 
terized by an exponentially decaying correlation function |6]]. For the specific case of CDW 
systems, one expects an exponential decay of the correlations between static fluctuations in 
the phase of the CDW order parameter j7|,|8| . 

In this paper, we report the results of a detailed high resolution x-ray scattering study 
of the static phase correlations of the Qi CDW in Ta-doped NbSe3. The experiments 
clearly show that phase fluctuations with high spatial frequency components are suppressed; 
therefore, the destruction of LRO by the random field is suppressed on length scales less than 
f« 75 A. These results are explained using the standard Ginzburg-Landau phase Hamiltonian 
with a Gaussian random pinning field [[F]|§. This new length scale is then related directly 
to the amplitude coherence length. We conclude with a brief discussion of the implications 
for physical models of the pinning mechanism. 

The obvious way to study the interaction between a CDW and a quenched random 
distribution of impurity atoms is to dope crystals with impurities and study how the ordering 
varies with impurity concentration and type. Over the last decade a large number of such 



studies have been reported PJlQ[]. Most of these studies have focussed on NbSe3, which 
undergoes two independent Peierls transitions to CDW states at T Pl « 145 K and at Tp 2 « 
59 K. The most widely studied dopant has been Ta, which is isoelectronic with Nb. 

Experimentally @,[J, the structure of the Qi CDW in Ta-doped NbSe3 is quite well 
described by a Ginzburg-Landau field theory |)ll|| . The effective Hamiltonian which describes 



the phase behavior of the Qi CDW at low temperatures can be written as 0] 

?U = f d d * [^ H • V0(x)) 2 + #(x)0(x)} , (1) 

where ip is the amplitude of CDW order parameter, </>(x) is the phase of the CDW order 
parameter, £ is the amplitude coherence length, and h(x) is the quenched random field. 

Using Tij, and the assumption that /i(x) obeys Gaussian statistics, the phase-phase cor- 
relation function in three dimensions is |7]|| /e*^( Xl )~^( X2 )] ) ~ e -s( x i- x 2) ) where 



9(x) =L /l (2^4iwf 1 - cos < k ' x »- (2) 



and L is the system size. At large separations, /e'^ xi ) ^ X2 ^\ ~ e l Xl X2 ^, where the 
correlation length is given by I = 32 \f . 

In order to provide a framework for discussing the experimental results, we first con- 
sider the consequences of random fields on the measured profiles. In an x-ray scattering 
experiment on a CDW system at low temperatures, the static structure factor <S(q) is pro- 
portional to the spatial Fourier transform of the phase-phase correlation function || . Thus, 
at low temperatures <S(q) exhibits Lorentzian squared fluctuations centered about the CDW 
satellite reflection positions. 

The experimentally measured profile is the convolution of <S(q) with the resolution func- 
tion of the diffractometer, 1Z(C)- If the sample has nontrivial mosaic structure, that must 
also be convolved into the resolution function. The form of Tl(C) can be calculated from first 
principles |12|]. To obtain a simple analytic form of the resolution function, we assume (i) 
that the line shape of the rocking curves of the monochromator and analyzer are Gaussians, 
(ii) that the half-widths of these curves are determined by the Darwin width of Si(lll), 



(Hi) that the angular profile of the wiggler beam is Gaussian, and (iv) that the angular 
divergence of the wiggler beam is that same as that of a bend magnet at the critical wave 
length. The assumed Gaussian line shapes allow the integrals to be done analytically. The 
result is then transformed from angle space to reciprocal space, and fit to an anisotropic 
Gaussian line shape. The half widths are W\\ = 2.08 x lO" 4 !" 1 and W ± = 2.06 x lO^A -1 . 
The long axis is rotated by 7.48° from the longitudinal direction. The ellipse describing the 
half height of R(C) is plotted in the inset to Fig. |l|. We will show below that the intrinsic 
scattering is, in fact, much narrower than the resolution of the diffractometer in both the 
q\\ (longitudinal) and the q z (out of scattering plane) directions. Therefore, the directions 
orthogonal to a transverse scan are integrated away and the dimensionality of the Fourier 
transform is effectively lowered from three to one . Thus, the observed line shape is given 
by the one- dimensional convolution of the intrinsic scattering with the resolution function. 
Since the one-dimensional Fourier transform of an exponentially decaying correlation func- 
tion is a simple Lorentzian, the observed scattering will be the convolution of a Gaussian 
resolution function and a Lorentzian structure factor. 

The experiments were carried out using a six-circle diffractometer on the X-25 wig- 
gler beam line at the National Synchrotron Light Source. The beam line optics include 
a toroidally focusing mirror, a double crystal Si(lll) monochromator, and a triple-bounce 
channel-cut Si(lll) analyzer crystal. Non- dispersive rocking curves of the analyzer produced 
Darwin-limited behavior. Our samples were single-crystal whiskers of NbSe3 lightly doped 
with Ta ( r(4k) — r fi ~ 45). The macroscopic sample dimensions were approximately 2 //m 
x 20 /iin x 5 mm. The crystal structure of NbSe3 is monoclinic and the real space lattice 
constants @ are a = 10.009A, b = 3.4805A, c = 15.629A, and (3 = 109.47°. We studied 



only the T Pl CDW whose wave vector is Qi = [0 Qi 0], where Qi ~ 0.243 and varies slightly 



with temperature |T^]. The samples were mounted using silver paint across a 3.5 mm hole 
in an alumina substrate and studied in transmission. Sample cooling was provided by a 
closed-cycle helium refrigerator. All of the data presented in this paper were collected at 
70 K. 



Figure [3] shows a scan in the b* direction through the [0 1 + Q\ 0] CDW satellite of 
sample #1. During a typical injection cycle, the stored electron beam current decayed from 
roughly 200 to 150 mA. The data was collected for a constant number of incident x-rays and 
is plotted as counts per sec assuming 100 mA in the storage ring. The half-width at half- 
maximum (HWHM) of the scan in the b* direction is 1.2 x 10~ 4 A _1 , which is 10% smaller 
than the value obtained from the cut through the a 'priori resolution function appropriate 
for a resolution-limited peak. This cut is indicated by the arrows in the inset to Fig. [I]. 
Thus, as claimed above, this scan is completely resolution limited. Figures ||] and ||| show 



scans through the [0 1 + Q\ 0] CDW satellite in the a* direction for two different samples. 
The Gaussian resolution function in the a* direction, which now includes the sample mosaic 
as measured at the [0 2 0] Bragg peak, is indicated by the dot-dashed line in the inset of each 
figure and the values of the HWHM, Wr, are listed in table |I[ The dashed lines in Figs. |2| 
and |3| are the best fit to the convolution of a simple Lorentzian and the Gaussian resolution 
for the two samples. The background level Ibg has been constrained to remain at or above 
the measured value of the background. A small asymmetry factor has been included. This 
factor improves the quality of the fit but does not affect the values of the other parameters. 
A casual inspection of Figs. |2] and ^ reveals two important points. First, the high brightness 
wiggler source associated with X-25 allowed us to obtain data over four decades of intensity, 
roughly two decades better than in previous measurements ||. Second, this extra dynamic 
range clearly reveals that the simple Lorentzian line shape overestimates the scattering in 
the wings. 

It is not surprising that the tails do not agree with a simple Lorentzian line shape. The 
exponential decay of the phase-phase correlation function is valid only at large separations, 
breaking down at small separations. The magnitude of this discrepancy can be made quan- 
titative by considering the Taylor series expansion of Eq. (|2|) about x = 0. 
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M-£t{j-z)* +0 W (3) 

In contrast to the asymptotic form, the full functional form of g(x) approaches the origin 



with zero slope. 

Previously, high resolution x-ray scattering studies of the structure of the Pt(OOl) sur- 
face [T7|] and high stage Br-intercalated graphite ]TS| have observed that it is necessary to 



approximate g(x) in a manner which preserves the functional form at both large and small 



distances |T!|. In our particular case, this can be done by letting g(x) = \/a 2 + x 2 /£. As 
before, the observed line shape is given by a one-dimensional convolution in reciprocal space 
of the static structure factor with the resolution function. The one-dimensional static struc- 
ture factor describing transverse scans through the CDW satellite produced by this model 
is p 

Va 2 + x 2 
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= KAa^- 2 + ql) (5) 

+ q± v J 

where K\(x) is the modified Bessel function of order 1. For small arguments, K\{x) ~ x _1 
and the simple Lorentzian form is recovered. For large arguments, K\{x) ~ -j= exp(— x) 
and the scattering decays exponentially. Qualitatively, introducing a cutoff at small x re- 
duces the scattering in the wings. 

The solid lines in Figs. |2| and |3] are the best fit to the convolution of Eq. (|5]) with the 
Gaussian resolution functions shown in the insets. The agreement between the fits and the 
data is striking. Quantitatively, for sample #1, the goodness-of-fit parameter \ 2 changes 
from 26.7 to 5.5 when the parameter a is allowed to rise up from zero. The best fit values 
for a and P. for the two samples studied are listed in Table Q. These best fit values of £ are 
consistent with those obtained by DiCarlo, et al. in Ref. |5|. 

As a test of our assertion that the intrinsic line shape is nominally a Lorentzian because 
we have integrated over the gn direction, we fit these data to the convolution of the resolution 
function with a Lorentzian raised to an arbitrary power rj. Similar to a cut off at small x, 
increasing the value of r] suppresses the scattering in the wings. However, the data is best 
described by a Lorentzian. Quantitatively, for sample #1, the best fit value of the exponent 



was T) « 1.09 ± 0.04 with \ 2 ~ 13, roughly twice that obtained using Eq. ([|), and excluding 
a Lorentzian to the | line shape. 

The principal result of our line shape analysis, the existence of a second length scale, is 
not sensitive to our assumption of a Gaussian resolution function. This is quite important 
since the true resolution function is certainly not expected to be Gaussian. Both dynamical 
diffraction [Q and thermal diffuse scattering [|16|] theory predict that a Bragg peak should 



have tails which decay only as q~ 2 . Due to the convolution, the observed scattering associated 
with a resolution function with wings would be enhanced relative to the Gaussian prediction. 
But, since the predicted effect is a decrease in the scattering, our essential result is insensitive 
to such tails. 

Our line shape analysis of the data demonstrates that high frequency phase fluctuations 
are suppressed in NbSe 3 . This suppression is consistent with the predictions of the Ginzburg- 
Landau field theory which specifically excludes phase fluctuations with spatial frequencies 
greater than l/£. The connection between our simple approximate form of g(x) and the 



exact formulation can be examined by equating the Taylor series expansion of \J a 2 + x 2 /£ 
with Eq. (|3p, yielding a = (37r/4)£. Thus, combining our two data sets, we find that the 
amplitude coherence length in the a* direction is £ a » = 34.4 ± 10. 3A. 

The observation that high frequency phase fluctuations are suppressed has implications 
for our physical picture of the pinning mechanism, supporting the following picture which 



is based on the notion of a Friedel oscillation [20|. The same divergence of the susceptibility 



which drives the Peierls transition causes the conduction electron density to respond to a 
point impurity by creating an oscillating charge density with wavenumber Ikp. Of course 
these oscillations cannot continue to infinite distance in a real sample. In general, we expect 
an exponential decay of the form e~ r ^ for these oscillations due to scattering of the conduc- 
tion electrons by other defects. The minimum value of the mean free path is, by definition, 
the amplitude coherence length £. Thus, a pinning site prefers not only a particular value of 
the phase, but that the phase be constant over a region of characteristic size £. In contrast 
to the physical description usually given for strong pinning in which the phase gradients are 



largely confined to a small region surrounding an individual impurity pl| , phase gradients 
are expelled from the neighborhood of impurity sites. 

Both strong and weak pinning theories are subsumed within the random field description; 
however, the structure of the CDW is determined primarily by the statistics of the random 
field, not its intensity. In our ongoing work studying the kinetics of CDW systems, the more 
general random field description is expected to be the most useful. 

The authors thank R.E. Thorne for providing the NbSe3 samples. This work was sup- 
ported by Cornell's Materials Science Center (NSF Grant No. DMR-88-1858-A02) and by 
the NSF (Grant No. DMR-92-57466). Additional support was provided by the AT&T Foun- 
dation. These data were collected on beam line X25 at the National Synchrotron Light 
Source, Brookhaven National Laboratory, which is supported by the U.S. Department of 
Energy, Division of Materials Sciences and Division of Chemical Sciences (Contract No. DE- 
AC02-76CH00016). 



REFERENCES 

[1] A.B. Harris, J. Phys. C7, 1671 (1974); G. Grinstein and A. Luther, Phys. Rev. B 13, 
1329 (1976). 

[2] Y. Imry and S.-K. Ma, Phys. Rev. Lett. 35, 1399 (1975). 

[3] L.D. Landau & E.M. Lifshitz, Statistical Physics, (Pergamon Press, New York) 1980. 

[4] E. Sweetland, C.-Y. Tsai, B.A. Wintner, J.D. Brock, and R.E. Thorne, Phys. Rev. Lett. 
63, 2841 (1989). 

[5] D. DiCarlo, R.E. Thorne, E. Sweetland, M. Sutton, and J.D. Brock, Phys. Rev. B (in 
press) . 

[6] P. Debye, H.R. Anderson, Jr., and H. Brumberger, J. App. Phys. 28, 679 (1957). 

[7] L.J. Sham and Bruce R. Patton, Phys. Rev. B 13, 3151 (1976). 

[8] K.B. Efetov and A.I. Larkin, Sov. Phys. JETP 45, 1236 (1978). 

[9] For comprehensive reviews of the electronic properties of CDW systems see P. Mon- 
ceau, in Electronic Properties of Quasi-One-Dimensional Materials (Reidel, 
Dordrecht, 1985), Pt. II, p. 139; G. Griiner, Rev. Mod. Phys. 60, 1129 (1988). 

[10] R.E. Thorne, Phys. Rev. B 45, 5804 (1992). See also, F. Levy and H. Berger, J. Cryst. 
Growth 61, 61 (1983). 

[11] H. Fukuyama and PA. Lee, Phys. Rev. B 17, 535 (1978); PA. Lee and T.M. Rice, 
Phys. Rev. B 19, 3970 (1979). 

[12] R.W. James, The Optical Principles of The Diffraction of X-rays, (Ox Bow Press, 
Woodbridge, 1962). 

[13] J.L. Hodeau, M. Marezio, C. Roucau, R. Ayroles, A. Meerschaut, J. Rouxel, and P. 
Monceau, J. Phys. Cll, 4117 (1978). 



[14] A.H. Moudden, J.D. Axe, P. Monceau, and F. Levy, Phys. Rev. Lett. 65, 223 (1990). 

[15] This integral is Eq. 8.338(2) in I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, 
Series, and Products, (Academic Press, New York, 1980), p. 938. 

[16] B.E. Warren, X-ray Diffraction, (Dover, New York, 1990). 

[17] D.L. Abernathy, S.G.J. Mochrie, D.M. Zehner, G. Grubel, and Doon Gibbs, Phys. Rev. 
Lett. 69, 941 (1992). 

[18] S.G.J. Mochrie, A.R. Kortan, R.J. Birgeneau, and P.M. Horn, Z. Phys. B 62, 79 (1985). 



[19] References |17| and pl| stressed that it was important to keep g(0) = 0; however, the 
present work clearly demonstrates that it is the functional form, not the absolute value, 
which must be maintained. An incorrect absolute value can simply rescaled away by a 
constant multiplicative factor. An incorrect functional form always produces an incorrect 
line shape. 

[20] Walter A. Harrison, Solid State Theory, (Dover, New York, 1980). 

[21] J.R. Tucker, W.G. Lyons, and G. Gammie, Phys. Rev. B 38, 1148 (1988); J.R. Tucker, 
Phys. Rev. B 40, 5447 (1989). 



10 



FIGURES 



FIG. 1. Scan in the b* direction through the [0 1 + Qi 0] CDW satellite of sample #1. 
Error bars represent counting statistics. Solid line is the best fit to a Gaussian. Inset shows ellipse 
depicting the half-height of the a priori resolution function. The arrows indicate the value predicted 
for the width. 



FIG. 2. Scan in the a* direction through the [0 1 + Qi 0] CDW satellite of sample #1. The 
dashed line is the best fit to the convolution of a Lorentzian with the Gaussian resolution. The 
solid line is the best fit to the convolution of Eq. (10) with the Gaussian resolution. The inset 
shows the central portion of the scan on a linear scale. 



FIG. 3. Scan in the a* direction through the [0 1 + Q\ 0] CDW satellite of sample #2. 
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TABLES 
TABLE I. Summary of parameters describing NbSe3 samples lightly doped with Ta (vr ~ 45). 

The errors on the parameters represent \-a confidence limits of the nonlinear least squares fit. 

Sample #1 Sample #2 

a 102 ± 30 A 60 ± 38 A 

4* 3 275 ± 175 A 3 050 ± 400 A 

W R (HWHM) 3.3 x 10~ 4 A" 1 1.32 x 10~ 3 A" 1 

Ibg 0.61 counts/sec 0.24 counts/sec 
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